Dynamical phase coexistence: 
a simple solution to the "savanna problem' 
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We introduce the concept of dynamical phase coexistence to provide a simple solution for a long- 
standing problem in theoretical ecology, the so-called "savanna problem". The challenge is to 
understand why in savanna ecosystems trees and grasses coexist in a robust way with large spatio- 
temporal variability. We propose a simple model, a variant of the Contact Process (CP), which 
includes two key extra features: varying external (environmental/rainfall) conditions and tree age. 
The system fluctuates locally between a woodland and a grassland phase, corresponding to the 
active and absorbing phases of the underlying pure contact process. This leads to a highly variable 
stable phase characterized by patches of the woodland and grassland phases coexisting dynamically. 
We show that the mean time to tree extinction under this model increases as a power-law of system 
size and can be of the order of 10,000,000 years in even moderately sized savannas. Finally, we 
demonstrate that while local interactions among trees may influence tree spatial distribution and 
the order of the transition between woodland and grassland phases, they do not affect dynamical 
coexistence. We expect dynamical coexistence to be relevant in other contexts in physics, biology 
or the social sciences. 
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I. INTRODUCTION 

Savannas are open systems that feature a continu- 
ous grass layer and a discontinuous tree layer. They 
appear across a wide range of climatic and ecologi- 
cal conditions, and are characterized by the stable, 
though variable, coexistence of two distinct ty pes of 
vegetation, trees and grasses Sarmiento 1984l |. This 



coexistence is dynamic in the sense that the den- 
sity of trees varies widely both in space and time, 
as recently confirmed by observations of cyclic tran- 
sitions between empty and dense tree distri butions 
' Moustakas et al. 20081 . IWieeand et al. 2006l| . Sa- 



vannas have been studied from both experimental 
and theoretical points of view, and have become an 
important subject of study in ecology. However, the 
origin, nature, and dynamics of savannas are not yet 
well understood. How is long-term coexistence of 
trees and grasses possible without the superior com- 
petitor taking over, as happens in other ecosystems 
(grasslands or woodlands)? This is one of the long 
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standing puzzles in theoretical ecology, commonly 
referred as the savanna problem. A tentative an- 
swer to this question is provided by niche models in 
which, assuming (soil, rain, etc) heterogeneities, each 
life form occupies the regions for which it is a su- 
perior competitor IWalter 197ll I Walker et al. 1981L 
IWalker and Nov-Meir 1982l |. This type of solu- 
tion is conceptually unsatisfactory and is, any- 
how, not supported by recent empirical observations 
iJeltsch et al. 2000l Iffiggins et al. 2000l iBond 20081 
IScholes and Archer 1997| . 
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van Wiik and Rodrig uez-Iturbe 20021 
Meyer et al. 2007 ]. For example, the storage effect 



hypothesizes that birth rate variability promotes 
species coexistence in communities of long-lived 
organisms, so that in an environment which is 
frequently adverse , a long life span buffers trees 
against extinction Chesson 2000j |. Other, similar in 
spirit, bu ffering me chanisms have also been proposed 
peltsch et al. 2000} . Following these studies, we 



2 



build a minimalistic savanna model that allows us to 
assess the contributions of the following features to 
long-term tr ee-grass coexistence : i) Variable weather 
conditions iJeltsch et al. 1996 H ipjgins et al. 20dol . 
Ivan Wiik and Rodriguez-Iturbe 200a |^ ii) Mean 
annual precipitation, which has been reported 
to enhance and lim it the maximum tree-cover 
Sankaran et al. 20051 ISankaran et al. 20081 . 

iBucini and Hanan 2007 1. iii) Adult tree longevity 
Higgins et al. 2000^ . and iv) Positive and 
negative local density-dep en dent tree inter- 
actions 



Jeltsch et al. 19961 



iMever et al. 20071 . 
Isee Calabrese et al. In pressf . Some exist- 

ing demographic savanna models, which in- 
clude these elements along with many oth- 
ers, can reproduce the ma in traits of real sa- 



vannas 
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Mever et al. 200^ However, these models do not 



clarify which ingredients are necessary to produce 
long-term coexistence and which are superfluous. 

Our goal is to construct a minimal stochastic model 
that explains phase coexistence in savannas. For that, 
we start from a cellular automata defined on a square 
grid, in which each site can be occupied either by one 
tree or grass, and that follows the dynamical rules 
of the standard contact process (occupied sites are 
trees and empty ones are grasses; see next section). 
Two phases characterize the system: grassland (grass 
only) and woodland (tree dominated). Dynamical 
phase coexistence appears when the system is driven 
by a random external driver (fluctuations in rainfall) 
which forces the system to visit the two phases ran- 
domly in time. In this way, coexistence appears in a 
broad region of parameter space without fine-tuning 
of parameters. We shall show that the range of dy- 
namical coexistence is much enhanced if the age of 
the trees is included, so that trees can typically en- 
dure harsh conditions. 

This dynamical coexistence is not indefinite as the 
rainfall fluctuations will eventually lead the system to 
the grassland or absorbing state, which is character- 
ized by the complete absence of trees. Unlike previous 
work on envir onmental fluctu ations and tree-grass 
coexistence [van Wiik and R odriguc z^turbe 20021 . 
IJeltsch et al. 199^ ~ Higgins et al. 200oj |^ we focus 
here on characterizing the timescales over which 
dynamical coexistence is robust. Our approach 
consists of studying stability by means of the mean 
life-time, that is, the mean time that takes the 
savanna to reach the final absorbing state. We show 
that this time diverges with system size as N" and 
can be enor mous for even moderately sized systems 



Leigh 1981 1 



Finally, let us remark that the generic phase co- 
existence observed in our model, i.e. coexistence 
occurring in a broad region of parameter space, is 



rarely find in other non-equilibrium model/systems 
and, in most case s, the domain of co existence is typ- 
ically small (see jMufioz et al. 2005j [ and references 
therein). Besides, given that the phenomenon of co- 
existence also happens in biological, physical and so- 
cial sciences, the mechanism proposed in this paper 
is expected to be useful in many other contexts. 



II. MODEL 

Consider the simple co ntact process (CP) 
[Harris 1974 iHinrichsen 2000l |. Each node {i,j) 
of a two-dimensional square lattice can be either 
occupied Zij — 1 (tree) or vacant Zij — (grass). 
The dynamics is as follows: a tree is randomly 
selected, and it is removed from the system with 
probability d, otherwise, with probability b, it 
generates an offspring, which is placed at a randomly 
chosen nearest neighbor (n.n) provided it was empty 
(i.e. short-range seed-dispersal). Every time a tree 
is selected, time t is increased by 1/N{t), where 
N{t) is the total number of trees in the system; t is 
increased in one unit, corresponding to one Monte 
Carlo (MC) step or "year", whenever all trees have 
been selected once on average. Fixing d, a phase 
transition appears at some critical value be- 

For b > be the system is in the active (the system 
has a non- vanishing density of trees and can dynam- 
ically evolve) phase - woodland -, while for 6 < 6c it 
is in the absorbing - grassland - phase with zero tree- 
density (as schematically illustrated in Fig. 1). The 
reason for this last name is that if the system reaches 
this absorbing (i.e. only-grass) state, the situation 
is irreversible, it remains indefinitely trapped in it. 
Since there is no spontaneous generation of trees, the 
dynamics ceases if there are only grasses. This type 
of transiti on is commonly te rmed an absorbing phase 
transition 



Hinrichsen 2000l| . 



Now, we introduce the following two extra ingredi- 
ents to model savannas: 

i) Fluctuating external conditions: We as- 
sume that the birth probability depends on exter- 
nal condi tions (mostly annual precipitation, but also 
fires, etc Sankaran et al. 2005j |). so b becomes a time- 
dependent random variable (see the zigzaging line in 
Figure 1). To account for possible temporal correla- 
tions in weather conditions we take a time-correlated 
(colored) noise as follows. With probability q a new 
value of b is extracted from an uniform random dis- 
tribution in [bmimbmax] at cach MC step, otherwise 
(with probability 1 — q) b is kept fixed. Here, we 
take bmin = 0, while bmax is the control parameter. 
This models weather cycles of typical length 1/q. In a 
similar way, one could include more periodic weather 
oscillations, as those induced by "El Nino", which 
lead to similar results. 

ii) Age: Field studies reveal that the mortality dis- 
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FIG. 1: Schematic representation of the system dynamics. 
The underlying pure model has two homogeneous phases: 
active (woodland) and absorbing (grassland), separated 
by a critical point (the transition could also be discon- 
tinuous as in Fig. |4|. The control parameter fluctuates 
in time and shifts from one phase to another (droughts 
and rainy periods). This, combined with long-living trees, 
prevents the system from reaching a homogeneous steady 
state. 



tribution of some savanna t rees is consistent acros s 
years and climatic regimes [Moustakas et al. 2006j |. 
and that they have long lifespans. To model this 
we define trees with an intrinsic age- variable, a{i,j) 
measured in years. The death probability is taken 
to be age-dependent: d d{a{i,j)). In particular, 
a random number, 77, is extracted from a Gaussian 
distribution of mean a„i and variance a (typically, 
ttm = 100 and a = 20). If a(i, j) > -q then the selected 
old tree is removed; otherwise nothing happens. 

Some other effects, such as density-dependent de- 
mographic rates, can also be easily implemented in 
the model. 



Density-dependence: Negative and positive lo- 
cal density-dependent death probabilities account for 
tree-tree competition and facilitation, respectively. 
Both of t hese effects have been repor ted to act in 
savannas 



see Calabrese et al. In Dressl |. To model 



competition between a young tree at site {i,j) (tree 
with age a(i,j) < Oest, where Oest is the establish- 
ment age) and its neighbors, we increase its death 
rate as a function of the number of its nearest neigh- 
bors trees NN(i,j) = J2(k i) where the sum is 
restricted to nearest neighbors sites of Thus, 
we take d{i,j) = 1 — exp [—NN{i,j)]. Contrarily, to 
model strong facilitation we consider d{i,j) — 1 for 
NN{i,j) = and 1, and d{i,j) = for NN{i,j) > 2, 
that is, trees born in sites with only a few occupied 
neighbors die with certainty, otherwise they survive. 



III. MODEL ANALYSIS AND RESULTS 

We first analyze the role of each ingredient we have 
added to the basic contact process separately. Taking 
a fluctuating b and a fixed death rate, i.e. switch- 
ing off" the age eff"ect, the system shifts randomly 
between the tendencies to be in the active (tree den- 
sity larger than zero) and in the absorbing (zero tree 
density) phase of the underlying pure model (notice 
the zigzagging trajectory in Fig. |T]). For b,nax > 
the system hovers around its critical point, while, if 
bmax < be (resp. bmin > be) the fluctuating system 
is in the absorbing (active) phase, i.e. it is a grass- 
land (woodland). In principle, if the time series of b 
happens to be adverse (i.e. b < be) for a sufficiently 
long time interval, any finite system falls into the ab- 
sorbing state; i.e. the system has variability but little 
resilience to long adverse periods. 

The effect of weather correlations is as follows: for 
< g << 1 the birth rate is constant for long peri- 
ods and the system typically jumps, every 1/q years, 
from a pure CP homogeneous state to another one 
and, therefore, when b takes a value smaller than 
be it falls ineluctably into the absorbing state. For 
g « 1, the rate of variation of b is very fast and the 
system does not have the time required to relax to 
any pure CP steady state, and reaches instead an av- 
eraged density value (see solid curve in Fig. [5]). For 
realistic intermediate values (e.g. q — 0.03), the sys- 
tem exhibits much larger oscillations (see the dashed 
and the dotted l ine in Fig. [2t which re semble those 
in real savannas [Moustakas et al. 2008 1. 

On the other hand, taking a fixed birth rate, q = 0, 
i.e. switching off the weather variability the 
model becomes a CP with age. For this, if trees die 
at a fixed given maximum age, a„i (i.e. with a = 0), 
the density is known to exhibit damped oscillations in 
time of period 2am, and to conve rge asymptoti cally 
to a homogeneous stationary value Gerami 2002j (see 
dashed-dotted line of Fig. [5]). Our model (with a > 0) 
exhibits analogous, though more variable, damped os- 
cillations, and converges either to the absorbing state 
or to an active homogeneous state. Separating these 
two regimes there is a phase transition. In this case 
the system is resilient but has little variability. Obvi- 
ously resilience grows with the maximum age. 

As we demonstrate below, the full model, which 
includes both age and correlated variable rainfall, ex- 
hibits variability and resilience (see the dashed and 
dotted curves in Fig. [2]); the larger the maximum age, 
the larger the resilience. The system fiuctuates lo- 
cally between the absorbing and the active phases, 
but it is able to preserve "islands" of the unfavored 
phase in a "sea" of the dominant phase, as is required 
for generic phase coexistence. This is the basic mech- 
anism of dynamical phase coexistence. 

The effect of density-dependence is as follows. 
Computer simulations show that, the underlying pure 
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FIG. 2: Time evolution of the tree-density in different 
cases. All curves except the dotted one {am = 30, cr = 6) 
are for am ~ 100, a = 20. The dashed-dotted curve 
{q — 0, i.e. model without variability) exhibits damped 
oscillations. The solid line {q = 1) shows small variability, 
while the dashed and dotted curves are for q — 0.03 (in- 
termediate variability). Notice that the one with smaller 
maximum age ~ 30; dotted curve) does not survive 
to an adverse period, while its analogous for large age 
{am ~ 100; dashed curve) does. 
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FIG. 3: Phase diagram for the underlying pure model 
{q — 0) in the case of strong facilitation; observe the 
discontinuous transition and the hysteresis loop. 
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FIG. 4: Stationary tree-density p" vs maximum birth 
probability bmax, averaged over surviving realizations for 
the model with competition, q = 0.1, a™, = ICQ, aest = 1, 
and system sizes A'^ — 40^ (diamonds). A'' = 80^ (squares) 
and N = 320^ (circles). 



model {q = and no age), either in the absence of 
density dependence and in the case of competition 
or weak facilitation, exhibits a continuous absorbing 
phase transition (as schematized in Fig. [I}. Instead, 
for strong facilitation the underlying transition can 
be discontinuous with a broad hysteresis loop (see 
Fig. [3]), implying that around the transition the two 
dynamically coexisting phases are very different. 

Fig. U shows the stationary tree density as a func- 
tion bmax {q = 0.1, Qm = 100, Uest = 1) illustrat- 
ing the existence of an active phase and a continuous 
phase transition in the full model with competition. 
The active phase of the full model is the phase of co- 
existence: the two phases of the pure model (grass 
and wood) coexist. A very similar continuous transi- 
tion is obtained in the case of strong facilitation; even 
if the underlying pure-model transition is discontin- 
uous, once varying conditions are switched on, the 
transition between the absorbing and the tree-grass 
coexistence phases becomes continuous. 

For illustration purposes. Fig. [S] shows snapshots 
of such a phase for different parameters at different 
times. Panels (a) and (b) correspond to the case 
of competition, and both have the same parameter 
values but look quite different, illustrating the large 
spatio-temporal variability. Panel (c) shows that, in 
the case of facilitation more compact clusters are ob- 
served (as justified by the underlying discontinuous 
transition) . These snapshots are visually very simila r 
to pictures of real savannas [Moustakas et al. 2008j ; 
notice the presence of irregularly distributed tree 
clusters of different sizes and shapes. Panel (d) corre- 
sponds to competition, but with slightly different pa- 
rameters. We can conclude that density-dependence 
controls the shape of emerging clusters, but it is not 
an essential ingredient for dynamical coexistence. 

Let us now scrutinize the stability of the active 
phase, where there is dynamical coexistence: how 
long does it persist through periods of adverse rainfall 
conditions? 

In a single realization, independent of the value of 
bmax and provided that bmin is smaller than the crit- 
ical values of the pure model, the density of trees 
reaches the absorbing state p = 0, that consists of a 
totally empty lattice where no more trees are born, 
and therefore the evolution stops. In other words, 
trees go extinct in a time that depends on the specific 
realization, due to stochastic fluctuations, and the 
size of the system. We have computed the mean life- 
time, T, that is, the average over many realizations 
of the time required to reach tree extinction, for sev- 
eral values of bmax and different system sizes TV (see 
Fig. [SJ. The value b'^^^ ~ 0.1045 at which r grows 
as T ~ (In A^)-^-^® [straight dashed line in Fig. [SJa)], 
agrees within error-bars with the point where the sta- 
tionary density p" goes to zero in Fig.|4j thus we take 
it as the active-absorbing phase transition point. In 
Fig. [n^b) we observe that, remarkably, in the active 



FIG. 5: Snapshots of the system for different versions of 4 
the model with q = 0.01 and hmax = 1. (a) and (b) corre- 
spond to the same realization with flest = 4, competition, 

and two different times. The system fluctuates locally 10^" " 



and globally between small and large densities, as in real 

savannas, (c) Facilitation, a^st = 4. (d) Competition, j^q' j^q^ j^q^ jq4 

Oest = 1. ^ 2 

N = L 



phase (circles) r diverges with TV as a power law A''" 
(with possible log corrections), where the exponent 
a increases with hmax (see inset). That is, the mean 
life-time increases very rapidly with system size, and 
therefore one expects that real (large) savannas sur- 
vive for very long periods. 

As an illustration, consider a 1000 hectare savanna, 
which corresponds to a 632 x 632 square lattice (as- 
suming, say, neighboring lattice sites are separated 
by 5 meter distance), taking h^aax — 0.111 not very 
deep into the active phase, then our model predicts 
that its expected time to extinction is of the order of 
10 million years. This leads to the conclusion that 
the coexistence of grassland and woodland phases is 
stable for extremely long time periods. Furthermore, 
we find that the exponent a increases with the max- 
imum tree age (not shown), so that age strongly 
stabilizes the active phase and increases the average 
savanna mean life-time time. 



IV. EMPIRICAL EXAMPLE 

To demonstrate that this simple model is a reason- 
able descriptor of savanna dynamics, we fit it to a 
paleoecological time series of tree-grass pollen ratios, 
Rdata{t), and lake depth levels, 6{t), from the Cres- 
cent Island Crater site in central Kenya [Lamb et al. 
2003, Vershuren et al. 2000]. The long-term nature 
of this time series (> 1000 years) allows us to exam- 
ine a much wider range of density fluctuations than is 
possible with short-term, directly observed tree den- 



FIG. 6: (a) Savanna mean life-time r vs the logarithm 
of its size A'^ on a log-log scale, for different values of 
bmax below (diamonds, bmax = 0.100 — 0.104) and above 
(circles, hmax ~ 0.105 — 0.115) the transition point b^^x — 
0.1045. At 6JL„^ r grows as 

(lniV)3-68 (straight dashed 
line), (b) r vs A on a log-log scale, r diverges as N" for 
bmax > bmax (clrcles) , wlth the exponent a proportional 
to bmax (inset), indicating that the active phase is stable 
when N — >■ cxd. 



sity data. The cost is that these paleoecological data 
sets rely on proxy measures (e.g. fossil pollen for 
vegetation, fossil midge and diatom assemblages for 
lake depth) and are thus necessarily much less pre- 
cise than direct measures of tree density and rainfall. 
We digitized figures 3B and 3C in [Lamb et al. 2003] 
to obtain the pollen ratio and lake depth time se- 
ries, respectively, using the program Engauge V 4.1 
I http://digitizer.sourceforge.net/ . As these time se- 
ries were both reconstructed at the same site, over 
the same time interval, they give us rough proxy mea- 
sures for both precipitation patterns and tree density 
at this site going back over 1000 years (see Fig. [7]). 

We ran simulations of the model with a birth prob- 
ability given by the expression b{t) = "fS(t) (where 
7 is a constant), that is, assuming that the amount 
of rainfall, and therefore b(t), is proportional to the 
lake's depth. We included neither age nor density ef- 
fects. We also considered that the amounts of pollen 
from trees and grasses are proportional to the num- 
ber of trees pL^ (with L = 100) and the area cov- 
ered by grass (1 — p)L^, respectively, thus we took 
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FIG. 7: Ratio of tree to grass pollen in the Crescent Is- 
land Crater core NC93 in Kenya (dashed-dotted line) and 
from numerical simulations of the model (solid line), on 
a 100 X 100 square lattice, with age = 100. The time 
series of the birth probability in the simulations (dashed 
line) was taken to be proportional to the depth of the lake 
Naivasha, located next to the Island Crater, in Kenya. 



Rnum{t) = — p) as the tree-grass pollen ratio 

in the simulations, where /? is a constant. Then, for 
given values of 7 and /3, we averaged over 10 inde- 
pendent realizations of the dynamics, and calculated 
the sum of squared deviations between the model and 
the tree-grass ratio data over the 1000 years period 

Comparing the model to the paleoecological data 
from the Crescent Island Crater site, we found by nu- 
merical direct search that 7min = 0.018 and Pmin = 
0.005 are the values that minimized the sum of 
squared errors. As can be seen in Fig. [71 the time 
evolution of the tree-grass pollen ratio from numeri- 
cal simulations (solid line) is strongly biased by the 
temporal function h{t) (dashed line), and it has only 
a qualitative agreement with the evolution of the real 
ratio (dashed-dotted line) . We have checked that the 
age of trees in the model does not affect the evolu- 
tion of the ratio very much, and that single realiza- 
tions are typically very similar, given that fluctua- 
tions coming from the system size are much smaller 
than fluctuations induced by the variation of h{t). 
We also note in Fig. [7] that the real-data pollen ra- 
tio also roughly follows the variation of the lake's 
depth. These results demonstrate that this correla- 
tio n between rainfall and pollen ratio, a lready noted 
in iLamb et al. 2003L IVerschuren et al. 2000] . IS cor- 
rectly described by our very simple model. This cor- 
respondence is especially encouraging given that we 
have manipulated only two parameters here to fit the 
model. 



V. SUMMARY 

Here we have shown that a simple extension of the 
contact process is capable of providing a robust and 



general explanation for tree-grass coexistence in sa- 
vannas. Specifically, fluctuating rainfall levels force 
the system to oscillate between a tree-dominated 
woodland (active) phase and a grass-only (absorb- 
ing) phase; a behavior we call dynamical phase co- 
existence. Dynamical coexistence is not permanent, 
but, as we have shown, it is expected to last for geo- 
logically significant periods of time. In other words, 
the mechanism discussed here facilitates coexistence 
over periods of time easily long enough to span the 
gaps between major disturbance events, such as ice 
ages, that can create and destroy savannas. 

Strictly speaking, the only essential ingredient for 
dynamical phase coexistence is the presence of ex- 
ternally varying conditions. Said another way, dy- 
namical coexistence is independent of the lifespans of 
trees, the degree of correlations in weather conditions, 
and the presence and nature of density-dependent lo- 
cal interactions among trees. Adding age-dependent 
tree death to the model greatly enhances the 
stability and robustness of coexistence. Adding 
weather correlations of intermediate magnitude in- 
creases the variability in tree cover and can help 
the model mimic the wide variability in tree cover 
that has been observed am o ng sites with similar rain- 
fall ISankaran et al. 200,1 H ucini and Hanan 20071 
ISankaran et al. 2008j |. Finally, adding local tree- 
tree interactions can allow the model to repro- 
duce a range of tree spatial patterns that have 
been observed in real savannas, including dense 
thickets of trees and more widely spaced and 
open configurations [see Calabrese et al. In pressl . 
iMoustakas et al. 2008j |. Thus these additional fea- 
tures affect the nature of several observable features 
of savannas but are not the primary drivers of tree- 
grass coexistence. 

This last point merits further discussion. While 
several stochastic savanna models have included 
varying weather condition s and many other 
factors IJeltsch et al. 1996L IJeltsch et al. 1998 . 
Higgins et al. 2000l [Mever et al. 2007 



van Wiik and Rodriguez-Iturbe 200*21 ^ none has, 
to our knowledge, clearly identified the minimal 
conditions that facilitate tree-grass coexistence. Our 
strategy of starting from the well-studied contact 
process and extending the model in a stepwise 
fashion has allowed us to identify the contribution 
of each model component to tree-grass coexistence 
and to other features observed in real savannas. 
We can therefore unambiguously state that fiuc- 
tuating external conditions alone are sufficient to 
facilitate long-term (though not indefinite) tree- 
grass coexistence. A further advantage of our 
minimalistic approach is that it does not require 
parameter fine tuning to achieve coexistence, as has 
been the case in o t her stochastic sava nna models 
jJeltsch et al. 19961 iHiggins et al. 2000f . Finally, 
unlike previous studies, we have been able to clearly 
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demonstrate the time scales over which dynamical 
phase coexistence persists. 

Despite the simplicity of our model, we have shown 
that it is capable of describing qualitatively a 1000 
year paleoecological dataset when it is given a quan- 
tity proportional to rainfall as a driving input. The 
key result here is that the model is able to repro- 
duce the correlation between weather patterns and 
tree-grass ratios. The rough quantitative agreement 
between the model and data is not surprising given 
that these data, which are based on proxy measures 
instead of direct observation, are relatively crude 
and that only two parameters of our dynamic, non- 
linear model, were tuned to achieve the fit. We 
note that though empirical data are frequently used 
in savanna modeling studies, such studies typically 
only conside r single snapshots in time, or very short 



time series [ Jcltsch ct al. 1999, Higgi ns et al. 2000l . 
Ivan Wiik and Rodriguez- I turbc 2002] . 

While we have focused here on savannas, the con- 
cept of dynamical phase coexistence is much more 



general than that. Variants of the Contact Process 
are used in many different fields including: Physics, 
Ecology, Epidemiology, Sociology, etc. Furthermore, 
it seems likely that many of the systems to which 
these models are applied may feature some degree 
of forcing by externally varying conditions. Thus, 
we expect that the concept of dynamical phase co- 
existence will find broad applicability in a range of 
scientific disciplines. 
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